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We study the performance of simple error correcting and error avoiding quantum codes together 
with their concatenation for correlated noise models. Specifically, we consider two error models: 
i) a bit-flip (phase-flip) noisy Markovian memory channel (model I); ii) a memory channel deflned 
as a memory degree dependent linear combination of memoryless channels with Kraus decomposi- 
tions expressed solely in terms of tensor products of Jf-Pauli (Z-Pauli) operators (model II). The 
performance of both the three-qubit bit flip (phase flip) and the error avoiding codes suitable for 
the considered error models is quantified in terms of the entanglement fidelity. We explicitly show 
that while none of the two codes is effective in the extreme limit when the other is, the three-qubit 
bit fiip (phase flip) code still works for high enough correlations in the errors, whereas the error 
avoiding code does not work for small correlations. Finally, we consider the concatenation of such 
codes for both error models and show that it is particularly advantageous for model II in the regime 
of partial correlations. 

PACS numbers: quantum error correction (03.67.Pp); decoherence (03.65. Yz). 

I. INTRODUCTION 

It is well known that one of the most important obstacles in quantum information processing is decoherence. It 
causes a quantum computer to lose its quantum properties destroying its performance advantages over a classical com- 
puter. The unavoidable interaction between the open quantum processor and its environment corrupts the information 
stored in the system and causes errors that may lead to wrong outputs. 

There are different methods for preserving quantum coherence. One possible technique exploits the redundancy in 
encoding information. This scheme is known as "quantum error correcting codes" (QECCs). For a comprehensive 
introduction to QECCs, we refer to [1]. Within such scheme, information is encoded in linear subspaces (codes) of 
the total Hilbert space in such a way that errors induced by the interaction with the environment can be detected 
and corrected. The QECC approach may be interpreted as an active stabilization of a quantum state in which, by 
monitoring the system and conditionally carrying on suitable operations, one prevents the loss of information. Another 
possible approach is the so-called "decoherence free subspaces" (DFSs) (also known as error avoiding or noiseless 
codes). For a comprehensive introduction to DFSs, we refer to [2]. It turns out that for specific open quantum 
systems (noise models in which all qubits can be considered symmetrically coupled with the same environment), it is 
possible to design states that are hardly corrupted rather than states that can be easily corrected (in this sense, DFSs 
are complementary to QECCs). In other words, it is possible to encode information in linear subspaces such that the 
dynamics restricted to such subspaces is unitary. This implies that no information is loss and quantum coherence is 
maintained. DFS is an example of passive stabilization of quantum information. 

The formal mathematical description of the qubit-environment interaction is usually given in terms of quantum 
channels. Quantum error correction is usually developed under the assumption of i.i.d. (identically and independently 
distributed) errors. These error models are characterized by memoryless communication channels A such that n- 
channel uses is given by A(") = A^". In such cases of complete independent decoherence, qubits interact with 
their own environments which do not interact with each other. However, in actual physical situations, qubits do 
interact with a common environment which unavoidably introduces correlations in the noise. For instance, there are 
situations where qubits in a ion trap set-up are collectively coupled to their vibrational modes 3J. In other situations, 
different qubits in a quantum dot design are coupled to the same lattice, thus interacting with a common thermal 
bath of phonons |4] . The exchange of bosons between qubits causes spatial and temporal correlations that violate the 
condition of error independence [5]. Memory effects introduce correlations among channel uses with the consequence 
that A^") ^ A®". Recent studies try to characterize the effect of correlations on the performance of QECCs [SHIP). 
It appears that correlations may have negative [7] or positive [S] impact on QECCs depending on the features of the 
error model being considered. 

Devising new good quantum codes for both independent and correlated noise models is a highly non trivial problem. 
However, it is usually possible to manipulate existing codes to construct new ones suitable for more general error 
models and with higher performances [111 I12| . Concatenation is perhaps one of the most successful quantum coding 
tricks employed to produce new codes from old ones. The concept of concatenated codes was first introduced in 
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classical error correcting schemes by Forney [T3]. Roughly speaking, concatenation is a method of combining two 
codes (an inner and an outer code) to form a larger code. In classical error correction, Forney has carried on extensive 
studies of concatenated codes, how to choose the inner and outer codes, and what error probability threshold values 
can be achieved. The first applications of concatenated codes in quantum error correction appear in |141 115| . In the 
quantum setting, concatenated codes play a key role in fault tolerant quantum computation and in constructing good 
degenerate quantum error correcting codes. For instance, Shor's degenerate code can be constructed by concatenating 
the three-qubit bit flip and phase flip codes. 

It is known that DFSs are efficient under conditions in which each qubit couples to the same environment (collective 
decoherence) while ordinary QECCs are designed to be efficient when each individual qubit couples to a different 
environment (independent decoherence) [16'. While none of the two error correction schemes is effective in the 
extreme limit when the other is, QECCs will still work for correlated errors [T71 [TS], whereas DFSs will not work in 
the independent error case Therefore, it would be of interest to study the performance of concatenated codes 
obtained by combining DFSs and QECCs for explicit noise models in the presence of partially correlated errors. 
Indeed, this will be one of the main purposes in our work. 

In this article, we study the performance of simple error correcting and error avoiding quantum codes together 
with their concatenation for correlated noise models. Specifically, we consider two error models: i) a bit-flip (phase- 
flip) noisy Markovian memory channel (model I); ii) a memory channel deflned as a memory degree dependent linear 
combination of memoryless channels with Kraus decompositions expressed solely in terms of tensor products of X-Pauli 
(Z-Pauli) operators (model II) . The performance of both the three-qubit bit flip (phase flip) and the error avoiding 
codes suitable for the considered error models is quantified in terms of the entanglement fidelity. We explicitly show 
that while none of the two codes is effective in the extreme limit when the other is, the three-qubit bit flip (phase 
flip) code still works for high enough correlations in the errors, whereas the error avoiding code does not work for 
small correlations. Finally, we consider the concatenation of such codes for both error models and show that it is 
particularly advantageous for model II in the regime of partial correlations. 

The layout of this article is as follows. In Section II, we briefly discuss about the concatenation technique in 
quantum coding, the DFSs and, the entanglement fidelity as a performance measure of quantum error correction 
schemes. In Section III, we evaluate the performances of the three-qubit bit flip code and suitable DFSs for error 
models I and II. The performance of the concatenated code for both models appears in Section IV. Finally, in Section 
V, we present our concluding remarks. 

II. ON CONCATENATION, DECOHERENCE FREE SUBSPACES AND, ENTANGLEMENT FIDELITY 

In this Section, we briefly discuss about the concatenation technique in quantum coding, the DFSs and, the 
entanglement fidelity. 

A. Quantum Concatenation Trick 

For the sake of simplicity we only use two layers of concatenation and consider single qubit encoding. The general- 
ization to arbitrary I layers of concatenation and multi-qubit encoding is relatively straightforward [20j . Assume the 
inner code (first layer) is a [ni, fci, di] stabilizer code Ci with generators Gi = {gj : i = I,..., ni — fci}, and the outer 
code (second layer) is a [n2, 1, ^2] stabilizer code C2 with generators G2 = {gj : j = f,.-., nj — l}. The concatenated 

dcf 

code C = Ci 0C2 is to map fci qubits into n — nin2 qubits, with code construction parsing the n qubits into ni blocks 
B (b) (6 = 1,..., ni) each containing n2 qubits. In other words, given a codeword |cinnor) for the inner code Ci, 

ICinnor) = ^1 '"^^^ ' 

where are basis vectors for Ci, the concatenated code C is constructed as follows. For any codeword |coutcr) for 

the outer code C2, 



(2) 
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with |zi...?„2) = •Si.-.'Si \in2)i replace each basis vector \ij) with ij = 0,..., ki — 1 for j = 1,..., n2 by a basis vector 
I ) in Ci , that is 

|Cconcatenatcd) ^ "n,..., i„2 l<^Jl) ® •■• ® I'^J-.a) • (3) 

ii,..., 

Further details on the construction of the stabilizer generators of C can be found in ^20] . As a final remark, we simply 
point out that the above mentioned construction produces a [nin2, fci, d] code with d > did2- As an illustrative 
example, consider a concatenated code that is based on I layers of concatenation of the same code, for instance the 
seven-qubit [7, 1, 3] CSS code. Here, an unencoded qubit is encoded into a block of seven qubits. Next each qubit is 
itself encoded into a block of seven qubits. Repeating this process I times produces the concatenated code in which 
one logical qubit is recursively encoded into the state of 7' physical qubits. Thus, such concatenated code is a [7', 1, d] 

code with d > 3'. Recall that a distance d QECC can correct t errors, where t = ^^^-^-^ ^rid [x] is the integer part of 

X. Since d > 3', we see that the number of errors that a concatenated code can correct grows exponentially with the 
number of layers of concatenation I. Therefore, concatenation can be a good coding trick as long as the error model 
at each encoding level has the same form, that is, the same Kraus error operators with possibly different amplitude. 



B. Decoherence Free Subspaces 



Following [2], we mention few relevant properties of DFSs. Consider the dynamics of a closed system composed of 
a quantum system Q coupled to a bath B. The unitary evolution of the closed system is described by the combined 
system-bath Hamiltonian -fftotj 

Htot = Hq®Ib + Hb®Iq+ Hint, Hint = ^ ® (4) 

a 

The operator Hq (i/g) is the system (bath) Hamiltonian, Iq (/g) is the identity operator of the system (bath), are 
the error generators acting solely on Q while Ba act on the bath. The last term in Q is the interaction Hamiltonian. 

A subspace "Hdfs of the total system Hilbert space "H is a decoherence free subspace (unitary evolution in "Hdfs 
for all possible bath states) if and only if: 

1. Ea IV") = Cq IV") with Cq G C, for all states \ip) spanning Hufs, a-nd for every error operator Ea in Hint- In other 
words, all basis states spanning Hops are degenerate eigenstates of all the error generators Ea', 

2. Q and B are initially decoupled; 

3. Hq \tp) has no overlap with states in the subspace orthogonal to Hdfs- 

To establish a direct link between QECCs and DFSs, it is more convenient to present an alternative formulation 
of DFSs in terms of the Kraus operator sum representation. Within such description, the evolution of the system Q 
density matrix is written as, 

PQ it) = TrB [U {pQ ® pb) U^] = ^^PQ (0) 4. (5) 

k 

where U — 6^^^*°** is the unitary evolution operator for the system-bath closed system and the Kraus operator A^. 
(satisfying the normalization condition) are given by, 

Ak = Vn {m\U\n) with ^ = ^Q^ (6) 

k 

where k = (n, m), |m) and \n) are bath states. It turns out that a iVDPS-dimensional subspace Hbfs of is a DFS 
if and only if all Kraus operators have an identical unitary representation (in the basis where the first A^dfs states 
span "Hdfs) upon restriction to it, up to a multiplicative constant. 
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where = ^Jn {m\Uc\n) and Uc = e with He = -ffg + iJinf Furthermore, Ak is an arbitrary matrix that acts on 

^DFS (with H — Hbfs ffi^ijFs) ^^-"^ cause decoherence there; J7g^^^^ is Uq restricted to 'Hdfs- Now recall that 
in ordinary QECCs, it is possible to correct the errors induced by a given set of Kraus operators {Ak} if and only if, 



Arfc/c 

■fe 



RrAk= [ Q ^ ^ , ),Vrandfc, (8) 



or, equivalently. 



where {Rr\ are the recovery operators. The first block in the RHS of ([s]) acts on the code space C while the matrices 
Brk act on where H = C (B ■ From ([t]) and ([s]), it follows that DFS can be viewed as a special class of QECCs, 
where upon restriction to the code space C, all recovery operators Rr are proportional to the inverse of the system Q 
evolution operator, 

i?.oc(<^«))\ (10) 
Assuming that (10 1 holds, from ([7| and ^ it also turns out that, 

AfeCxC/g'^'\ (11) 

upon restriction to C. Furthermore, from (?]) and (|9|, it follows that ^tk' — 9*k9k' ■ However, while in the QECCs case 
7fcfe' is in general a full-rank matrix (non-degenerate code), in the DFSs case this matrix has rank 1. In conclusion, a 
DFS can be viewed as a special type of QECC, namely a completely degenerate quantum error correcting code where 
upon restriction to the code subspace all recovery operators are proportional to the inverse of the system evolution 
operator. As a side remark, in view of this last observation we point out that it is not unreasonable to quantify the 
performance of both active and passive QEC schemes by means of the same performance measure. 



C. Entanglement Fidelity 

We recall the concept of entanglement fidelity as a useful performance measure of the efficiency of quantum error 
correcting codes. Entanglement fidelity is a quantity that keeps track of how well the state and entanglement of a 
subsystem of a larger system are stored, without requiring the knowledge of the complete state or dynamics of the 
larger system. More precisely, the entanglement fidelity is defined for a mixed state p = '^iPiPi =tr^^ lip) {ip\ in 
terms of a purification Itp) G H Hr to a reference system Hr. The purification lip) encodes all of the information 
in p. Entanglement fidelity is a measure of how well the channel A preserves the entanglement of the state "H with its 
reference system T-Ir. The entanglement fidelity is defined as follows |21j . 

-F(p, A)'^=^'(^|(A®/«J(|V')(^|)|V'), (12) 

where \ip) is any purification of p, is the identity map on M {'Hr) and A(8)/Kij evolution operator extended 

to the space H (8) "Hr, space on which p has been purified. If the quantum operation A is written in terms of its Kraus 
operator elements {Ak} as, A (p) = J^k ^kpA^, then it can be shown that 

T (p, A) = ^ tr (Akp) tr (4p) = ^ |tr {pAk)f . (13) 

k fe 

This expression for the entanglement fidelity is very useful for explicit calculations. Finally, assuming that 

A : M (H) 3 p ^ Aip) =J2 AkpAl eMiU), dimcH (14) 

fe 

and choosing a purification described by a maximally entangled unit vector G H H for the mixed state p — 
dhkn^-H ' obtain 



T ( A 1 = 



fe 



^Em^i'- (15) 



A in (14) 



The expression in ( 15 ) represents the entanglement fidelity when no error correction is performed on the noisy channel 
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III. THREE-QUBIT BIT FLIP CODE AND DFS 

In this Section, we consider two correlated error models. Error correction is performed by means of the three-qubit 
bit flip code and a suitable decoherence free subspace. The code performance measure used is the entanglement 
fidelity. 

Model I. The first model is a bit flip noisy quantum Markovian memory channel A^")(p) (model I). In explicit terms, 
we consider n qubits and Markovian correlated errors in a bit flip quantum channel. 



A(»)(p) J (g) ... (g) 



'A.J 



Zl,..., in—0 



where Aq Ai X are Pauli operators. Furthermore the conditional probabilities Pi^.\i^ are given by, 



with, 



X! P^„\^„-lPi,^-l\^^-2■■■Pi2l^^P^l =1- 
ii,..., in—0 



(16) 



(17) 



(18) 



To simplify our notation, we may choose to omit the symbol of tensor product "(g" in the future, Ai^C 
Ai^...Ai-^. Furthermore, we may choose to omit the bar "|" in Pii^\i^ and simply write the conditional probabilities as 

Pikij ■ 

Model II. The second quantum communication channel A|"^ (p) (model II) that we consider is a memory quantum 
channel defined in terms of a linear combination of simple memoryless channels with Kraus decompositions expressed 
in terms of bit-flip error operators. The coefficients of such combination are dependent on the memory parameter /z. 
In explicit terms, we consider the following channel A^""* (p). 



A(") (p)'^=ifao(/i)A(") {p)+a, (^) A^") (p) , 



(19) 



flip channel). Finally, A^"'' (p) describes n-uses of a memoryless channel whose Kraus decomposition is characterized 

n i 

only by weight-0 and weight-n error operators with amplitudes (1 — p) ^ and [1 — (1 — p) ] ^ , respectively, 

A^") (p) =^ (1 -p)" >+ [1 - (1 -p)"]Xix2...X"-iX"pXix2...X"~iX". (20) 

In this Section, QEC is performed via the three-qubit bit flip code and a suitable decoherence free subspace. Although 
the error models considered are not truly quantum in nature, from this preliminary work we hope to gain useful 
insights for extending error correction techniques to quantum error models in the presence of partial correlations. The 
performance of quantum error correcting codes is quantified by means of the entanglement fidelity as function of the 
error probability p and degree of memory fi. 



where oq (/i) =^ 1 — /i, ai (fi) =^ /i, Ag"' (p) = A'"-* (p) in ( 16 1 in the limiting case of /i = (n-uses of a memoryless bit 



A. The three-qubit bit flip code 

1. Model I 



Error Operators. In the simplest example, we consider the limiting case of ( 16 1 with n = 3, 

1 1 

A'^Hp) = P^3\i2Pi2\^lP^l A,^A,^A,^pA\^A\^A\^ , with ^ Pi^\^2P^2\^lP^l = ^- (21) 

il, i2, 13=0 il, 12, 13=0 



the following error operators, 

A ^ {A^,.., A,} with A(3)(p) '^j^ApA^l and, j^^A = h 



Substituting (|T7| in (pi]), it follows that the error superoperator A associated to channel (|2T| is defined in terms of 

8x8- (22) 



fc=0 



A;=0 



6 



In an explicit way, the error operators {A'q,.., A'-^} are given by, 

A'o = \/p^I^ «) (8) /^ A[ = \/^X^ (g) /2 (g) /3, A'2 = \/p^I^ (g) (g) /3, 



'/^g)X^ 

= ^p|i^^/ig)X2g)X3, = ^pf^Xig)X2g)X3, (23) 

~(3) 

where the coefficients p^, for k — 1,.., 7 are given by, 

~(3) 2 ~(3) ~(3) ~(3) 

Po = PooPo, Pi = PooPwPo, P2 = PoiPwPo, P3 = PooPoiPi, 

~(3) ~(3) -(3) ~(3) 2 fnA\ 

Pi = PioPilPo, Ph = PoiPwPi, Pe = PoiPiiPi, P) = PiiPi, (24) 

with, 

Po = (1 - P) , Pi = P, Poo = ((1 - (1 - P) + m) , 

Poi = (1 - A*) (1 -P) , Pio = (1 - /^)p, Pii = ((1 - Ai)p + m) ■ (25) 

Encoding Operator. Consider a three-qubit bit flip code that encodes 1 logical qubit into 3-physical qubits. The 
codewords are given by, 

|0) |0) ® |00) = 1000) |0l) , |1) |1) ® |00) = |100) ""--^'^ |110) |iii) . (26) 

The operator f/cNOT ^-'^^ CNOT gate from qubit i to j defined as, 

f^c^OT = \ W + ^0 ® I' + - ^0 ® ■ (27) 
Finally, the encoding operator [/one such that C/onc |000) = |000) and [/cnc |100) = |111) is defined as, 

C/cnc (f/cNOT ® ° (f^CNOT ® • (28) 

Correctable Errors and Recovery Operators. The set of error operators satisfying the detectability condi- 
tion 23 , PcA'^Pc = ^A'^Pc, where Pq = |0_l)(Ol|+ [II) is the projector operator on the code subspace 
C = Span{\QL) , is given by, 

ctcctablc ~ {A'„,A[,A'„A':„A'^,A',,A',}CA. (29) 

The only non-detectable error is A'j. Furthermore, since all the detectable errors are invertible, the set of correctable 
errors is such that ^Jorj-octabic-^corroctabio is detectable. It follows then that, 

•^correctable — {^0' ^1' ^2; ^3} — -'^detectable 

C A. (30) 
The action of the correctable error operators ^correctable on the codewords |0l) and |1l) is given by, 

|0l) ^ A'„ \0l) = VpFiOOO) , A[ \0l) = ^|100) , A', \0l) = Vp^|010) , A', \0l) = V^FlOOl) 

lu) ^ = VpF 1111), .4; iu) = ^ 1011), A^iu) = ^1101), a:,Iu) = VpF 1110). (31) 

The two four-dimensional orthogonal subspaces and V^^ of "Hj generated by the action of ^correctable on |0l) 
and |1l) are given by, 

V"^ = S-panjlw"^) = |000) , = |100) , = |010) , jw^^) = |001) }, (32) 
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and, 



V^^ S'pan{|t;J^) = |111) , = |011) , jug^) = |101) , jw^^) = |110)} , (33) 

respectively. Notice that © V^^ — li^. The recovery superoperator H. ^ {Ri} with I — 1,..,4 is defined as [18], 



(34) 



i=0 



where the unitary operator VJ is such that V/ \vl^) = jii) for « S {0, 1}. Substituting (32) and (33l into (34|, it 
follows that the four recovery operators {Ri, i?2, -R3, R4} are given by, 

Rl = |0L)(0L| + |lL)(lL|,i?2 = |0L)(100| + |lL)(011|, 



i?3 = |0L)(010| + |lL)(101|,i?4-|0L)(001| + |lL)(110| 



(35) 



Using simple algebra, it turns out that the 8x8 matrix representation [i?;] with / = l,--,4 of the recovery operators 
is given by, 



[Rl] = Ell + Ess, [R2] — E12 + Es7, [R3] — Ei3 + Esq, [Ra] — E14 + Ess, 



(36) 



where Eij is the 8x8 matrix where the only non-vanishing element is the one located in the ij-position and it equals 
1. It follows that TZ O {Ri} is indeed a trace preserving quantum operation since, 



R^Ri = hxs- 



(37) 



The action of this recovery operation TZ on the map A'-'^' (p) in ( 22 ) leads to. 
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A(l_ (p) ^ (7^oA(3)) (p) ''^^ ^ ^ (i?, A^) p (R,A',)^ , 



(38) 



fc=o 1=1 



Entanglement Fidelity. We want to describe the action of T^.oA^'^) restricted to the code subspace C. Therefore, we 
compute the 2x2 matrix representation [i?;^^ of each RiA'i^ with I = 1,.., 4 and k = 0,.., 7 where. 



dcf / (Oi|i?,A'jOi) (Oi|ii;,A;,|ii) 



(39) 



Substituting (311 and (35) into (39 1, it turns out that the only matrices with non- vanishing trace are given 

by. 



[RiA'^]\c = ypo 



A3) ^ 1 



1 



1 



[fi3^2]|C ^ \/P2 

Therefore, the entanglement fidehty T^-J {fi, p) defined as 

^(f,) (M, P) (1/2x2, 7^oA(3)) ^ ^ EE {[RiA',] 



;J|C 



fe=0 1=1 



results. 



~(3) , ^(3) . --(3) . -.(3) 



(40) 



(41) 



(42) 



The expression for J-^.^^ (/x, p) in (41) represents the entanglement fidelity quantifying the performance of the error 
correction scheme provided by the three-qubit bit flip code here considered. The quantum operation TJ.oA^'^) appearing 
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in (41) is defined in equation (38) and the recovery operators Ri are explicitly given in (35). The action of RiA'/, in 



(41 



is restricted to the code space C defined in ( 26 ) 



Substituting ([24|) and ([25| into (|42j) , we finally obtain 

j:-^^) p) = ^2 ^2p3 _ ^ p) + ^ (_4p3 ^ Qp2 _ + (2p3 _ 3p2 ^_ 1) (i^odel I). (43) 
Notice that for a vanishing degree of memory the entanglement fidelity becomes, 

•^bft (0, P) = 2p3 - 3p2 + 1. (44) 

Remarks on the coding for phase flip memory channels. The code for the phase flip channel has the same character- 
istics as the code for the bit flip channel. These two channels are unitarily equivalent since there is a unitary operator, 
the Hadamard gate H, such that the action of one channel is the same as the other, provided the first channel is 
preceded by H and followed by W j24] . 

^phasc dcf ^ ^bit ^ ^ (1 _ p) p + pZpZ, (45) 

where, 

A*^'* (p) {l-p)p + pXpX and, H ^ \ \ y (46) 

These operations may be trivially incorporated into the encoding and error-correction operations. The encoding for 
the phase fiip channel is performed in two steps: i) first, we encode in three qubits exactly as for the bit flip channel; 
ii) second, we apply a Hadamard gate to each qubit, 

|0) 1000) 1000) ^ |0l) !+ + +>, |1) |100) |111) ^ |1l) = I ) , (47) 



where. 



1 



%/2 VI -1/ Vli 



The unitary encoding operator J/Pjf"^ is given by Ui^^'''' =^ H'^^ o C/g^^*,with U^^l defined in hA. Furthermore, in 



(48) 



the phase flip code, the recovery operation is the Hadamard conjugated recovery operation from the bit flip code, 

j^phaso dof ^(g)3 jjbit^(gi3 

2. Model II 



We consider A^"-* (p) in ( 19 ) with n = 3. Technical details will be omitted. They can be obtained by following the 
line of reasoning presented for the three-qubit bit flip error correction scheme applied to the model I. It turns out 
that, 

(^, p) = p (-3p^ -f' 6p2 - 3p) + {2p^ - 3p2 + 1) (model II). (49) 

B. DFS for bit flip noise 

1. Model I 



Error Operators. Consider the limiting case of ( 16 ) with n = 2 qubits and correlated errors in a bit flip quantum 
channel, 

1 

Ai^^p)"^' E P^,\^,P^AA^,^A,Jp{A,^A,)\ (50) 
il, 12=0 
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The error superoperator A associated to channel (50) is defined in terms of the following error operators, 

A ^ {A'„,.., A',} with A(2)(p) '='^^',K^ and, ^ A^U;, = h^,. (51) 

k=0 fc=0 

In an explicit way, the error operators {Ag,.., A'^} are given by, 

'^ri^ra _ , ^ 7-2 4' _ , /^(2) rl ^ ^-2 /l' _ , /^(2) ^ v2 



^0 = VPqI^ ® = yPiX^ ® ^2 = yP2l^ ® ^3 = yn 'X^ ® (52) 

where the coefficients p^'^-' for fc — 0,.., 3 are given by, 

-(2) -(2) -(2) -(2) ,ro\ 

Po = PooPo, Pi = PwPo, P2 = PoiPi, P3 = PiiPi- (53) 
Encoding Operator. We encode our logical qubit with a simple decoherence free subspace of two qubits given by [2] , 

|0) \0l) = and, |1) = \-+) , (54) 



with |±) defined in (48). As a side remark, we point out that a suitable DFS encoding for a phase flip Markovian 
correlated noise model is given by |0) — > |0l) = |01) and, |1) — > = |10). 

Correctable Errors and Recovery Operators. The set of error operators satisfying the detectability condi- 
tion 23|, PcA'i^Pc = ^A'^Pc, where Pq = |0_l)(Ol|+ is the projector operator on the code subspace 
C = S'pan {|0l) , |1l)} is given by, 

-4dctcctablc = {^0' ^3} ^ (55) 

Furthermore, since all the detectable errors are invertible, the set of correctable errors is such that 
-^Lrrcctabic-^corroctabiG is detectable. It follows then that, 

-^correctable •^detectable 

c A. (56) 
The action of the correctable error operators -4correctabie on the codewords |0l) and |1l) is given by, 

|0i) ^ A', = 4W\+-) , ^3 |0l) = -\/pF|+-) , |1l) -> A', \U) = , A', lU) = -Vpf . 

(57) 

The two one-dimensional orthogonal subspaces V"^ and V^^ of 'H\ generated by the action of y^correctabic on |0l) 
and are given by, 

V"^ = S-panlji;""-) = 1+-)} and, V^^ fi-pan { jw^^^ ) = |-+) } . (58) 

Notice that V^^ 7^ 'H2- This means that the trace preserving recovery superoperator TZ is defined in terms of 

1 

one standard recovery operator Ri and by the projector R±^ onto the orthogonal complement of V*^, i. e. the 

1=0 

part of the Hilbert space T-L\ which is not reached by acting on the code C with the correctable error operators. In 
the case under consideration, 

2 

def I 



^1 = l + -> + . fi± = E l''^) ' (59) 

5-1 

where Iks)} is an orthonormal basis for (V°^ ® V^^)^. A suitable basis B, „ ^ nj- is given by. 

Therefore, TZ O Ri_} is indeed a trace preserving quantum operation, 

R\Ri + R\R^ = h^i. (61) 
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The action of this recovery operation 72. with R2 = R±_ on the map A*^^^ (p) in (51 1 yields, 



3 2 

A(,^_ (p) ^ (72oA(2)) (p) '^^f ^ ^ [RiA'^) p iR,A',)^ , (62) 



Entanglement Fidelity. We want to describe the action of TZoK^"^^ restricted to the code subspace C. Therefore, we 
compute the 2x2 matrix representation [/ijAy |^ of each RiA'^. with I — 1, 2 and k — 0,.., 3 where, 



. d_of { {OlIRiA'M {Ql\RiA',\1l) \ 




by, 



FIG. 1: Threhold curves for model I: concatenated code (dashed line), DFS (thin solid line). 
Substituting (52) and (59 1 into (63), it turns out that the only matrices with non- vanishing trace are given 



[^1^0] 



I 
^" I 1 



3J|C 



^3 I 1 



(64) 



(2) 

Therefore, the entanglement fidehty J'j^pg (/i, p) defined as, 



results, 



(2) 



DFS 



(65) 



(66) 



The expression for J'^pg (/i, p) in (65) represents the entanglement fidelity quantifying the performance of the error 
correction scheme provided by the error avoiding code here considered. The quantum operation T^-oA^^-* appearing in 
(65) is defined in equation (62 1 and the recovery operators Ri are explicitly given in (59). The action of RiA'^. in (65) 



is restricted to the code space C defined in (54) 



Substituting ( 25 ) and ( 53 1 into ( |66[ ) , we finally obtain 

J'Sls (m, P)=l^ (-2p2 + 2p) + (2p2 - 2p + 1) (model I). 



(67) 



We point out that error correction schemes improve the transmission accuracy only if the failure probability V (/i, p) 
is strictly less than the error probability p [20] , 
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dcf 



p) - 1--F(a., p)<p. 



(68) 



In view of (68), we can determine threshold curves p.{p) that allow to select the two-dimensional parametric region 
where error correction schemes are useful. For instance, considering the model I in absence of correlations, it follows 
that the three-qubit code is effective only ii p < 0.5. However for such values of the error probability, the DFS 
considered does not work for fi approaching zero. The threshold curve for the DFS for the model I appear in Figure 
1. The DFS works only in the parametric region above the thin solid line in Figure 1, while the three-qubit bit-flip 
code works for all values of the memory degree /U when the error probability is less than 0.5. 



2. Model II 



We consider A^"' (p) in (19) with n = 2. Technical details will be omitted. They can be obtained by following the 



line of reasoning presented when studying the error avoiding code applied to the model I. It turns out that, 

J'Sls (a*: p) = (- V + 2p) + (2p2 ~2p+l) (model II). 



(69) 



Notice that the entanglement fidelities in (67) and (69) are equal. Finally, the threshold curves for the three-qubit 



bit flip and DFS for the model II appear in Figure 2. The DFS works only in the parametric region above the dashed 
line while the three-qubit code works only below the thin solid line. 




FIG. 2: Threshold curves for model II: DFS (dashed line), three-qubit bit flip code (thin solid line), concatenated code (thick 
solid Une). 



IV. CONCATENATED CODES 



A. Model I 



Error Operators. Consider the limiting case of ( 16 ) with n — 6 qubits and correlated errors in a bit flip quantum 
channel, 



(70) 



ii,..., 16=0 



The error superoperator A associated to channel ( 70 ) is defined in terms of the following error operators, 

A ^ {A'o,.., A',,} with A(«)(p) ""^^ ^'kP^k and, ^ ^A'^ - /64x64. (71) 



fe=0 



fe=0 
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The error operators in the Kraus decomposition ( 71 1 are 2° = 64 

6 



E(!)-2^ (72) 



k=0 

where (^) is the cardinality of weight-fc error operators. 

Encoding Operator. We encode our logical qubit with a concatenated subspace obtained by combining the deco- 
herence free subspace in (54) (inner code, Cdfs — Cjnncr ) with the three-qubit bit flip code in (26) (outer code, 
Cbit = Couter)- We obtain that the codewords of the concatenated code C = Cofs ° Cbit are given by, 

|0l) = ^ (1000000) - 1000111) + 1111000) - 1111111)), 

|li) = i (1000000) + 1000111) - 1111000) - 1111111)). (73) 

As a side remark, we point out that a suitable concatenated code for the case of a phase flip Markovian correlated 

noise model is given by |0) — > |0l) = |+ + + ) and, |1) — |1l) = | h ++). 

Correctable Errors and Recovery Operators. Recall that the detectability condition is given by PcA'/^Pq = A^j^Pc 
where the projector operator on the code space C is Pc — |0l) (Ol| + |1l) Observe that, 

PcA'^Pc = |0i) (Oil + (IlI + {Il\A'M |1l) (Ol| + |1l) (1l| • (74) 

Therefore, it turns out that for detectable error operators we must have, 

{Ol\A'^\Ql) = {IM^l) and, (Oi|A',|U) = {IM'M = 0. (75) 

In the case under consideration, it follows that the only error operators (omitting for the sake of simplicity the proper 
error amplitudes) not fulfilling the above conditions are proportional to, 

X^X^X^ and X'^X'^X^. (76) 

For such operators, we get 

{Ql\X^X''X^\Ql) = 1, (IlI^'^'^'IU) = -1, and, {{)l\X^X'> X''\Ol) = -1, (1l|^^^'X'|1l) = 1- (77) 

Therefore X^X^X^ and X'^X^X^ are not detectable. Thus, the cardinality of the set of detectable errors -Adotectabio is 
62. Furthermore, recall that the set of correctable errors -Acorrectabie is such that A^^^^^^^^^y^^^Acomctahie is detectable (in 
the hypothesis of invertible error operators). Therefore, after some reasoning, we conclude that the set of correctable 
errors is composed by 32 error operators. The correctable weight-O, 1 and 2 correctable error operators are (omitting 
the proper error amplitudes), 

{-^iweight-o ' ^"^^ ^^ ^^ ^^iwoigM-i > c^^) 

and, 

[X^X"^, X^X"^, X^X^, X'^X^, X'^X^, X^X^, X^X^, X^X^, ^^^^}„^ight-2 ' (''9) 
respectively. The correctable weight-4 errors are, 

X^X^X-^X'^ x^x^x^x^, x^x^x^x^, x^x^x^x^, x^x^x-^x^^ 
X^X^X^'X^, x^x^x^x^, x^x'^x^x^, x^x'^x^x^ 

Finally, weight-5 and weight-6 error operators arc given by. 



(80) 



weight -4 



{X^X^X^X'^X'\ X^X^X^X'^X^, x^x'^x^x^x^, x^x^x'^x-'x^, x^x^x^x'^x^, x'^x^x^x^x^ I . , , 

I ' ' ' ' ' J wcignt-5 ' 

(81) 

and, 

{X^X'^X^X^X'^X'^\ (82) 

I J wcight-6 ' ^ ^ 



13 



respectively. The action of the correctable errors on the codewords in ( 73 1 is such that the Hilbert space H2 can be 
decomposed in two 32-dimensional orthogonal subspaces V"^ and V^^ . In other words, "Hf = "^^^^ © where 



= Span { \vl-^,) = ^=A', \0l) } and, V'^ = Span { - ^=A', |U; 



(83) 



with A'l. e y^corroctabic Vfc ~ 0,..., 31 (numbering the correctable error operators from to 31). Notice that (^vl^ \vl'/ ^ — 
Skk'Sij, with k, k' £ {0,..., 31} and i, j e {0, 1} since. 



»L ; , \JL 

, M6) Mfi) 

^ yPk yPk' ' 



^L\K^k'\3L 



-^'kk'^ij — ^kk'Sij, 



(84) 



where we have used the fact that the square (Hermitian) matrix aj^j,, equals y P'^^pf) 5kk' ■ The recovery superoperator 
n O {Ri} with I = 1,.., 32 is defined as |18], 



where the unitary operator Vi is such that Vi \v\'") — jii) for i £ {0, 1}. Notice that, 



(85) 



(86) 



If turns out that the 32 recovery operators are given by, 

^^+1 - = (|0l) (Oil + |1l) (IlI) 

with ^ G {0,..., 31}. Notice that TZ O {Ri} is a trace preserving quantum operation since. 



(87) 



32 



32 



E^?^' = E Ei^-)(-M Eb- 



1=1 



1 = 1 \ IL 



32 



= E E l«r)(^^b-^)(-rhEEI^^)(^/"h^' 



64x64) 



;=1 iL 



Since with I = 1,..., 32 and e {0, 1} is an orthonormal basis for Hf- Finally, the action of this recovery 

operation TZ on the map A^^^ (p) in ( 71 ) leads to, 



2"-! 32 



(89) 



fc=0 i=l 



Entanglement Fidelity. We want to describe the action of TZoK^^^ restricted to the code subspace C. Recalling that 



A[ — A'l , it turns out that. 



{iL\Ri+i^k\jL) = ^ (*l|Ol) (oiiAfA'jji) + (iiiu) 



(90) 



We now need to compute the 2x2 matrix representation [_R;Ay|^ of each RiA'^ with / ~ 0,.., 31 and k — 0,.., 2^ — 1 
where, 



[Ri+iAk\\c 

For I, k — 0,.., 31, we note that [i?;+i74y|g becomes, 

, [Ol\A[^A'^\Ol 
[Ri+iA'k] • ' 



,1 del f {Ol\Ri+iA'^\Ol) {Ql\Ri+iA',\Il) 
Il\Ri+iA'^\Ol} {1l\Ri+iA',\1l) 



lL\AWk\lL 



PlOlk 



1 
1 



(91) 



(92) 
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while for any pair {I, k) with I = 0,.., 31 and k > 31, it follows that, 

{Ol\Ri+iA'^\Ol) + {Il\Ri+iK\Il) = 0. (93) 
We conclude that the only matrices with non-vanishing trace are given by with I = 0,.., 31 where. 



Therefore, the entanglement fidelity J^;^onc (/^, p) defined as 



becomes. 



2°-l 32 



:f!:± (/., p) 7^oA(«)) - ^ E E ([^'^^] 



fe=0 ( = 1 



(94) 



(95) 



1.00 



0^5 - 



0^ ■- 



0.55 ■- 



0.50 




FIG. 3: Entanglement fidelity vs. memory parameter /i with p = 10 ^ for model II: concatenated code (dashed line), DFS 
(thin solid line) and three-qubit bit flip code (thick solid line). 



•^conc (m, p) = PaoPo + 2pooPioPQ + pIoPoiPio (4po + Pi) + "^ploPaiPioPa + PmPoiPioPiiPo + 3pooPoiPioPo 



+4^PoiPwPiiPo + ^PoiPioPiiPi + PooPoiPioPliPi +PoiPiopIi [po + ^Pi) + '^PwPtiPo + pliPi- (96) 



Substituting (|25| into ( 96 ) , we finally get 



T^L (P, P) = 



24p^ - 24/ 



(40p^ - 120/ + 130/ - 60/ + 10/) + 



{-80p^ + 240/ - 264p'' + 128/ - 26/ + 2p) - 
(80/ - 240/ + 252/ - 104/ + lOp^ + 2p) + 



~H (-40/ + 120p'^ - 112/ + 24/ + 12p2 - Ap) 



3/ - 24p^ + 18/ + 4/ - Qp^ + 1) . (97) 



The threshold curves for the concatenated code defined in (73) for the model I appear in Figure 1. It turns out that 



the concatenated code does not work in the region delimited by the two dashed lines. We emphasize that in view of 
equations (43), (67) and (97), it turns out that the concatenated code outperforms the bit-flip code in regions with 



very high memory parameter values and outperforms the DFS in regions with low memory parameter values. In 
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particular, the relevance of the concatenation trick shines where the weaknesses of the inner and outer codes (yet not 
concatenated) are more pronounced, that is in regions with both high error probability and low memory parameter 
values. It is within this area that we can identify the region where the concatenated code outperforms both the inner 
and the outer codes. However, there is a big parametric region characterized by intermediate values of the degree of 
memory (see Figure I) where the concatenation trick docs not work well for model I. On the contrary, we will show 
that the quantum coding trick is particularly useful for the model II in the presence of partial correlations. 



B. Model II 



We consider A^"'' (p) in (19) with n = 6. Technical details will be omitted. They can be obtained by following the 



line of reasoning presented when studying the concatenated code applied to the model I. It turns out that. 



•^conc(M, P) 



2Ap^ - 18/ - + 6p2) + (8p^ - 24p^ 



1 



(98) 



The threshold curve for the concatenated code defined in (73) for the model II appears in Figure 2. It follows that 



while none of the two codes is effective in the extreme limit when the other is, the three-qubit bit flip (phase flip) code 
still works for correlated errors, whereas the error avoiding code does not work in the absence of correlations. The 
concatenated code works everywhere except below the thick solid curve in Figure 2. Here there is a parametric region 
characterized by partial correlations and delimited by the dashed and thin solid lines where only the concatenated code 
is effective. Furthermore, from ( [49| , (|69| and (|98|, it follows that the concatenated code is especially advantageous 
for the model II for partially correlated error operators. For the sake of clarity, in Figure 3 we plot the entanglement 
fidelities (69) (thin solid line), (98) (dashed line) and (49) (tick solid line) for p = 10~^. For such value of the error 
probability, the error avoiding code only works for /i > 0.44 (threshold value obtained from Figure 2), the three-qubit 
bit flip code only works for /i < 0.30 (threshold value obtained from Figure 2) while the concatenated code with 
entanglement fidelity given in (98) is is efficient for any value of /i G [0, 1]. 



V. FINAL REMARKS 



In this article, we studied the performance of simple error correcting and error avoiding quantum codes together with 
their concatenation for correlated noise models. For model I, a bit-flip (phase-flip) noisy Markovian memory channel, 
we have applied both the three-qubit bit flip (phase flip) and a suitable error avoiding code. The performance of the 
codes was quantified in terms of the entanglement fidelities in ( 43 ) and ( 67 ) . The performance of the concatenated 
code applied to model I appears in (97). In Figure 1, we have plotted the parametric regions where error correction 
is effective. We have presented a similar analysis for the model II, a memory channel defined as a memory degree 
dependent linear combination of memoryless channels with Kraus decompositions expressed solely in terms of tensor 
products of X-Pauli (Z-Pauli) operators (model II). The performance of the codes was quantified in terms of the 
entanglement fidelities in (49) and (69). The performance of the concatenated code applied to model II appears in 
(98). In Figure 2, we have plotted the parametric regions where error correction schemes are effective. 

Our analysis explicitly shows that while none of the two codes is effective in the extreme limit when the other is, 
the three-qubit bit flip (phase flip) code still works for correlated errors, whereas the error avoiding code does not 
work in the absence of correlations. Finally, our flnal finding leads to conclude that the concatenated code in ( 73 ) is 
particularly advantageous for model II in the regime of partial correlations (see Figure 3). 
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